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Abstract 

Isoscaling, where ratios of isotopes emitted from two reactions exhibit an 
exponential dependence on the neutron and proton number of the isotope, 
has been observed over a variety of reactions including evaporation, strongly 
damped binary collision, and multifragmentation. The conditions for isoscal- 
ing to occur as well as the conditions when isoscaling fails are investigated. 
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With the availability of rare isotope beams as well as detection systems that can resolve 
the masses and charges of the detected particles, isotope yields become an important observ- 
able for studying nuclear collisions of heavy ions [1,2]. This additional freedom on isospin 
asymmetry allows one to study the properties of bulk nuclear matter that are affected by 
the nucleon composition of the nuclei such as the isospin dependence of the liquid gas phase 
transition of nuclear matter [3,4] and the asymmetry term [5-8] in the nuclear equation of 
state. To minimize undesirable complications stemming from the sequential decays of pri- 
mary unstable fragments, it has been proposed that isospin effects can best be studied by 
comparing the same observables in two similar reactions that differ mainly in isospin asym- 
metry [4,6,8]. If two reactions, 1 and 2, have the same temperature but different isospin 
asymmetry, for example, the ratios of a specific isotope yield with neutron and proton num- 
ber N and Z, obtained from system 2 and system 1 have been observed to exhibit isoscaling 
i.e. exponential dependence of the form [4,5]: 

R 21 (N, Z) = Y 2 (N, Z)/Y X {N, Z) = C ■ exp(N -a + Z-(3) (1) 

where a and (3 are the scaling parameters and C is an overall normalization constant. 
We adopt the convention that the neutron and proton composition of reaction 2 to be more 
neutron-rich than that of reaction 1. 

Figure 1 illustrates the isoscaling property observed with the fragments produced in the 
central collisions of 124 Sn + Sn and 112 jf>n+ 112 Sn reactions [4]. The N and Z dependence 
of Eq. p] becomes most apparent if R 2 i (N, Z) is plotted versus N or Z on a semi- log plot as 
shown in the left panels. Isotopes of the same elements are plotted with the same symbols. 
Odd-Z nuclei are represented by open symbols while the even-Z nuclei are represented by 
closed symbols. In the upper left panel, the isotopes for each element Z appear to lie on 
one line and the resulting slopes would then be a. The dashed (for odd-Z elements) and 
solid (even-Z elements) lines are best fits to the data points with one common a values for 
all the elements. In this case, a = 0.361. Similarly, plotting R 2 i(N,Z) against Z for all 
isotones would provide a common slope /3 for each N. This is demonstrated in the lower 
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left panel of Figure 1. The best fits (dashed and solid lines) yield a value of (3 — —0.417. 
The positive (negative) slopes of the lines in the top (bottom) panel arise from the fact that 
more neutron-rich (proton-rich) nuclei are produced in the more neutron-rich (proton-rich) 
system which represent the values in the numerators (denominators) of the ratios in Eq. 1. 

Alternatively, the data in the left panels can be displayed compactly as a function of one 
variable, either N or Z, by removing the dependence of the other variable using the scaled 
isotope or isotone functions [6]: 



For the best fit value of —0.417), S(N) for all elements lies on a single line on a 
semi-log plot as a function of N as shown in the upper right panel of Figure 1. Alternatively, 
S(Z) of all isotones lies on a single line on a semi-log plot as a function of Z for the best fit 
value of a(— 0.361) as shown in the lower right panel. Both solid lines shown in the right 
panels have the same exponential dependence on N and Z as the corresponding lines in the 
left panels. The agreement between the data and the lines is excellent, verifying the scaling 
relation of Eq. 1-3. In general, the fit for S(N) is better than S(Z); this may reflect the 
influence of Coulomb interaction which may not be well approximated by Eq. 1 [6,8]. On 
the other hand, S(N) is affected mainly by the differences in neutron chemical potentials or 
the neutron separation energies. These latter factors may be governed by the differences in 
the symmetry energies in the two systems [8]. 

In a recent survey of heavy ion induced reactions, isoscaling appears to be manifested 
in a variety of nuclear reactions including deeply inelastic collisions, evaporation and multi- 
fragmentation over a wide range of incident energies [6]. In this article, we will perform a 
comprehensive exploration of many reactions and examine conditions under which isoscal- 
ing is observed and others where it is not. We will also demonstrate how isoscaling can be 
restored even when two systems have different temperatures. 



S(N) = R 21 (N,Z) xexp(-/3Z). 



(2) 




(3) 
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I. Deeply Inelastic Collisions 

In the 1970's, the Deeply Inelastic Collision (DIC) phenomenon was discovered when 
heavy ions were used to bombard heavy targets in an effort to create superheavy elements 
[9, 10]. Products from the DIC exhibit characteristics which can partly be attributed to 
compound nuclei decay and partly to multinucleon transfer reactions depending on the 
incident energy and detection angles. In ref. [6] the isotope yield ratios of 16 + 232 Th 
and 16 + 197 Au reactions at incident energy of 137 MeV and 9 = 40° have been found to 
exhibit isoscaling behavior. From the literature, we have selected four additional systems to 
illustrate the compliance or noncompliance of the scaling behavior in DIC . Each pair of the 
chosen reactions use the same projectile at the same incident energy and detect the isotopes 
at the same laboratory angles. The main differences are the targets. Figure 2 shows the 
relative isotope ratios, R 2 i(N,Z) for the four systems: a.) 16 + 232 Th [9] and 16 + l97 Pb 
[10] at incident energy of 137 MeV and 9 = 40° (upper left panel), b.) U N + 100 Mo and lA N 
+ 92 Mo at 97 MeV and 9 = 25° (upper right panel) [9], c.) 22 Ne + 232 Th and 22 Ne + 94 Zr 
at 173 MeV and 9 = 12° [9] (lower left panel), and d.) 16 + 232 Th and 16 + 197 Au at 315 
MeV and 9 = 40° (lower right panel) [10]. Isotopes of the same elements are plotted with 
the same symbols using the same convention as Figure 1, open circles, closed circles, open 
square, closed squares, and open diamonds for Z=3, 4, 5, 6, and 7, respectively. The solid 
and dotted lines connect isotopes of the same elements (solid lines for even Z element and 
dashed line for odd Z elements). Scaling similar to Equation 1 is observed for the isotope 
ratios plotted in the upper panels. When the product nuclei are detected at very forward 
angles, such as the isotope ratios from the 22 Ne induced reactions shown in the lower left 
panel, scaling is not observed. When the incident energy is raised to 20 MeV per nucleon, 
target dependence is much weaker than at lower energies [10] and production of isotopes 
at forward angles is consistent with fragmentation of the projectile and shows no target 
dependence. In that case, a = 0,(3 = and i?2i(iV, Z) 1 is observed (lower right panel). 

Figure 2 summarizes that isotopic scaling is reasonably well respected at low incident 
energies (E/A < lOMeV") and at angles backward of the grazing angles but poorly respected 
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at forward angles. The situation at higher energies is not clear. Isoscaling may have been 
observed with very small values of a and (3. The positive observation of isoscaling can be 
understood as follows: Backward of the grazing angles, it is often assumed that equilibrium 
is established between the orbiting projectile and target. In such cases, the isotopic yields 
follow the " Q gg -systematics" [9,10], in which the primary isotope yield of the projectile-like 
fragment depends mainly on the Q-value of the mass transfer and can be approximated by 

Y(N,Z) oc exp((M P + M T - ATp - M^)/T) (4) 

where Mp and My are the initial projectile and target masses, and M' p and M' T are the final 
masses of the projectile- and target-like fragment. Here, T has a natural interpretation as the 
temperature, but is not always assumed to be so. Applying charge and mass conservation, 
and expressing explicitly only the terms that depend on N and Z, one can write R 2 %(N, Z) 
as 

R 21 {N, Z) oc exp[(BE(N 2 -N,Z 2 -Z)- BE(N 1 -N,Z X - Z))/T], (5) 

where Zi and iVj are the total proton and neutron number of reaction i. BE is the bind- 
ing energy of a nucleus. Expanding the binding energies in Taylor series, one obtains an 
expression of the form 
BE(N 2 -N,Z 2 -Z)- BE(N l -N,Z 1 -Z)^a-Z + b- N + c-Z 2 + d-N 2 + e- ZN, 

(6) 

where a, b, c, d, and e are constants from the Taylor expansion. Evaluating Eq. |^ within the 
context of a liquid drop model, one finds that the second order terms are of the order (1/A), 
where A is the mass number, relative to the first order terms. The leading order parameters, 
a and b can be interpreted as the differences of the neutron and proton separation energies 
for the two compound systems, i.e. a = —As n and b = —As p . Equation |5| can then be 
approximated as 

R 2l (N, Z) oc exp[(-A Sn • N - As p ■ Z)/T\. (7) 
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Eq. |7| confirms the earlier studies which showed that the symmetry contribution in As n 
of the various isotopes associated with the same element, is approximately linear in the 
number of neutrons transferred. Similarly, As n shows a linear dependence on the charge 
transferred due to Coulomb-barrier effects [10, 11]. Comparison of Eqs. [I] and [7| reveals that 
the difference in the average separation energies, As n /T and As p /T , plays a corresponding 
role to the fitting parameters of a and (3. From the binding energy expansion in Eq. |], one 
expects that Eq. [7| becomes less accurate and eventually breaks down leading to a failure in 
isoscaling when the range of fragment masses considered becomes large. 

To explore how good is the approximation of using the nucleon separation energies, we 
calculate i? 2 i obtained with Eq. [7] as well as using the exact expression of Eq. |. We use the 
two parent systems that describe the 16 + 197 Au and 16 + 232 Th reactions corresponding 
to (Ni, Z\) = (126,87) and (N 2 , Z 2 ) = (150,98). The deviations, R 21 {Eq. $/R 2 i(Eq. § are 
plotted in Figure 3 for Li, Be, B, C, N and O isotopes as a function of the neutron excess (N- 
Z). Within each element, the deviations assume a parabola shape with the minima located 
at the neutron rich side of the N=Z line. Over the range of nuclei measured experimentally, 
\N — Z\ < 3, the overall deviations are less than ±5%. However, this comparison suggests 
that isoscaling will break down for isotopes with large Z and for nuclei with extreme isospin 
asymmetry. It appears likely that first order deviations from the scaling behavior can be 
corrected using functions similar to the parabola shown in Figure 3. 

II. Compound nucleus decay 

The scaling behavior for fragments evaporated from an excited compound nucleus has 
been discussed in ref. [6]. The measured isotope ratios for 4 He + 116 Sn and A He + 124 Sn 
collisions at E/A = 50MeV [12] are plotted in Figure 4 using the same convention of Figure 
1 and 2. At back angles (8 = 160°, left panel), the isotope ratios of different elements have 
similar slopes and adjacent elements are separated nearly equi-distance from each other, 
typical behavior of isoscaling as evidenced by the best fit solid and dashed lines. However, 
at forward angles (8 = 12°) where contributions from pre-equilibrium processes become 
significant, isoscaling is not well respected as shown in the right panel. Not only do different 
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elements have different slopes but also the distances between adjacent isotones vary greatly 
[12]. In general, there is a tendency for the slope to become steeper as the fragment mass 
is increased, consistent with the heavier elements emitted at lower temperature. However 
for the carbon isotope yields, the trend is actually not monotonic with N, indicating a clear 
failing of isoscaling. More detailed discussions on the forward angle data can be found in 
Ref. [12] suggesting that the failure of isoscaling may arise from non-equilibrium emission. 

The origin of isoscaling for evaporation process follows similar derivations involving the 
expansion of the binding energies in the Taylor series as described in previous section, 
resulting in a formula similar to Eq. (|7|). 

R 21 (N, Z) ex exp[(-A s „ + A/*) • N - (As p + A/* + A$) • Z)/T\. (8) 

where /* and /* are the neutron and proton excited free energy and $ is the electrostatic 
potential at the surface of a nucleus. A full derivation of Eq. (§) can be found in Ref. [6]. 

III. MULTIFRAGMENTATION PROCESS 

The scaling phenomenon was first observed in multifragmentation process in the central 
124 Sn + 124 Sn and 112 Sn + 112 Sn collisions[4] as demonstrated in Figure 1 and discussed in 
detail in the introduction section. To obtain guidance of how the nuclei yield ratios may 
be systematized, we examine the dependence of the isotopic yields within the equilibrium 
limit of the Grand- Canonical Ensemble [13,14]. In this case predictions for the observed 
isotopic yield are governed by both the neutron and proton chemical potentials, \x n and 
fjL p and the temperature Tj, plus the individual binding energies, BE(N, Z), of the various 
isotopes [11,12]. 

Yi(N, Z, = Fi(N, Z, Ti) ■ exp(N • /^/T; + Z • // P /Ti) • exp(BE(N, Z)/T ; ) (9) 

The factor F^N, Z, Tj) includes information about the secondary decay from both particle 
stable and particle unstable states to the final ground state yields. If the two reactions have 
the same temperature, (Tj = T) the binding energy terms in Eq. (9) cancel out. If one 
further assumes that the influence of secondary decay on the yield of a specific isotope is 
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similar for the two reactions, i.e. F±(N, Z,T) ~ F2(N, Z,T), then we obtain an equation in 
the same form as Eq. (1): 

R 21 (N, Z) = C- exp(N • A/x n /T + Z • A^/T) (10) 

where a = Afi n /T and (3 = A/j, p /T reflect the differences between the neutron and proton 
chemical potentials for the two reactions and C is an overall normalization constant. Afi n 
and A/ip correspond to As n and As p of Eq. [7| Simulations adopting microcanonical and 
canonical [8] multifragmentation models show that Eq. [X0| is respected. Recent SMM model 
calculations [15] indicate that fj, n and s n are closely related (f/, n ~ — s n +/*) for < T < 
3MeV, where the decay configurations are mainly binary, but the connection between fi n and 
s n becomes increasingly weak as the role of multifragment decay configurations becomes 
important. These calculations also verify the insensitivities of isoscaling to the effect of 
sequential decays [8]. 
IV Mixed Systems 

The isoscaling described by Eq. [I] relies on the emission mechanisms for the fragments 
in each reaction being described statistically with some common effective temperature and 
that distortions from secondary decays cancel[4,6,8,16,17]. The exhibition of the systematic 
trends does not imply that both reacting systems proceed with the same reaction mecha- 
nism. This point was demonstrated in Ref. [18] where isotopic yields of fragments produced 
in central Au + Au multifragmentation process at E/A = 35MeV [19] can be related approx- 
imately via isoscaling to those produced in lower multiplicity evaporation process produced 
in Xe + Cu reactions at E/A = 30MeV[20]. Isoscaling arises because the temperatures for 
the two reactions are nearly the same i.e. T\ w T 2 > [21] even though the emission mechanisms 
in the two reactions differ significantly [19,20]. 

For reactions which differ mainly in temperatures, isoscaling is also destroyed because 
the binding energy terms in Eq. ^ do not cancel even if the effect of sequential decays can 
be neglected. 

R 21 (N, Z) = C ■ exp(N • a' + Z ■ ft) exp(BE/T 2 - BE/Ti) (11) 
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where a'=a — k{i n2 and fl'=(3 — kfi p2 ; k = 1/Ti — 1/T 2 . While the new scaling parameters 
a' and f3' are related to a and (3, they do not have simple physics interpretations. The 
left panel of Figure 5 shows the R 2 \ ratios extracted from isotope yields of l24 Sn + 124 Sn 
and 4 if e + 124 Sn. Even though the comparison of the two Sn + Sn reactions and the two 
alpha induced reactions exhibit isoscaling as seen in Figs. 1 and 4, respectively, there is no 
observable scaling in these systems with different temperatures. 

Isoscaling could be restored if R 2 \ (N, Z) in Eq. |IT] is multiplied by the Boltzmann factor 
with binding energy and temperatures, exp(k * BE(N, Z)). Previous studies suggest that 
the temperature of the multifragmentation reaction of 12A Sn + 124 Sn collision is about 5 
MeV [22] and the temperature of the evaporation reaction of 4 He + 124 Sn is about 3 MeV 
[12], we obtain k pa 0.12. In Figure 4, R 2 i(N, Z) exp(0.12*BE) obtained from the same data 
plotted in the left panel exhibit very nice systematic behavior shown in the right panel. The 
restored isoscaling is clearly demonstrated by the dashed and solid lines which are the best 
fits through the data points with a' = 0.939. 

Currently, most temperature measurements depend on selected isotope yields e.g. the 
T iso (HeLi) depends on the yields of 3 ' 4 He and 6 ' 7 Li [2,23,24] and T iso (CLi) relies on the yields 
of n ' 12 C and 6 ' 7 Li [23]. Discrepancies in temperature measurements have been observed be- 
tween Tj<j (HeLi) and Tj SO (CLi) [23]. Furthermore, temperatures derived from excited states 
(T ex ) disagree with isotope yield temperatures (Tj so ) obtained from central collisions at inci- 
dent energy greater than 35 MeV [23,24]. Such discrepancies could arise if the light charged 
particles with Z < 2 are emitted early and/or the emitting sources are not thermalized 
[20]. With the temperature corrected isoscaling (Eq. |ll|), the internal consistency of the 
temperature measurements and the degree of thermalization as a function of incident energy 
can be investigated further using all the isotopes measured instead of selected isotopes. 

In summary, we have surveyed many reactions with different reaction mechanisms. We 
found that isoscaling occurs if both reactions can be described by statistical reaction mecha- 
nisms and that the temperatures of both reactions are nearly the same. However, isoscaling 
does not yield any information about the reaction mechanisms. In order to draw correct 
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conclusions from isotopic measurements, it is therefore absolutely essential to obtain ad- 
ditional experimental information that elucidates the underlying reaction mechanism. If 
the temperatures for both reactions are different, isoscaling can be restored with appropri- 
ate temperature corrections. This work was supported by the National Science Foundation 
under Grant Nos. PHY-95-28844 and PHY-96-05140. 
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Figure Caption: 

Figure INuclei yield ratios are plotted as a function of N (top panels) Z (bottom 
panels) for central 12A Sn + 124 Sn and ll2 Sn + 112 Sn collisions at E/A = 50MeV. The lines 
in the upper left panels correspond to best fits of different elements with one common slope. 
Similarly, in the bottom left panels, the lines correspond to fits of the same isotones. In the 
top right panel, the scaled isotopic ratio, S(N) (Eq. 2) is constructed using ft = —0.417. 
Similarly, in the bottom right panel, the scaled isotone ratio, S(Z) defined in Eq. 3, is 
plotted function of Z using a = 0.361. 

Figure 2 Relative isotope ratios for four systems: a.) 16 + 232 Th [9] and 16 O + 197 Pb 
[10] at incident energy of 137 MeV and 9 = 40° (upper left panel), b.) U N + 100 Mo and 
« jV + 92 Mo at 97 MeV and 9 = 25° [9] (upper right panel), c.) 22 Ne + 232 Th and 22 Ne + 
94 Zr at 173 MeV and 9 = 12° [9] (lower left panel) and d.) 16 + 232 Th and 16 + 197 Au 
at 315 MeV and 9 = 40° [10] (lower right panel). 

Figure 3 Deviations in approximating R 2 \ calculated with Eq. [7| from values calculated 
with Eq. [5] plotted as a function of the neutron excess for lithium to oxygen isotopes for the 
deeply inelastic reactions of [N u Zi) = (126,87) and (N 2 ,Z 2 ) = (150,98). 

Figure 4 Relative isotope (Z=3-6) ratios for A He + 116 Sn and A He + 124 Sn systems 
emitted at backward (left panel) and forward angles (right panel). See Figures 1 and 2 for 
symbol conventions. 

Figure 5 Left panel: Disappearance of isoscaling in reactions with different temperature. 
Right panel: Isoscaling is restored if the binding energy terms in the isotope ratio is corrected 
for the temperature difference. See Eq. ([□]). 
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